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Yerevan Physics Institute, 375036 Yerevan, Armenia 
(Dated: June 2, 2015) 

We present the final results for the deuteron spin structure functions obtained from the full data set 
collected with Jefferson Lab's CLAS in 2000-2001. Polarized electrons with energies of 1.6, 2.5, 4.2 and 
5.8 GeV were scattered from deuteron ('^^NDs) targets, dynamically polarized along the beam direction, 
and detected with CLAS. From the measured double spin asymmetry, the virtual photon absorption 
asymmetry Af and the polarized structure function gf were extracted over a wide kinematic range 
(0.05 GeV^ < < 5 GeV^ and 0.9 GeV < W < 3 GeV). We use an unfolding procedure and a 

parametrization of the corresponding proton results to extract from these data the polarized structure 
functions and g" of the (bound) neutron, which are so far unknown in the resonance region, W <2 
GeV. We compare our final results, including several moments of the deuteron and neutron spin structure 
functions, with various theoretical models and expectations as well as parametrizations of the world data. 
The unprecedented precision and dense kinematic coverage of these data can aid in future extractions of 
polarized parton distributions, tests of perturbative QCD predictions for the quark polarization at large x, 
a better understanding of quark-hadron duality, and more precise values for higher-twist matrix elements 
in the framework of the Operator Product Expansion. 

PACS numbers: 13.60.Hb, 13.88.+e , 14.20.Dh 
Keywords: Spin structure functions, nucleon structure 


I. INTRODUCTION 

One of the enduring goals in the field of hadron physics is 
a complete picture of how the fundamental particles of the 
standard model, quarks and gluons, make up the structure 
and the properties of the nucleon. Among other observ¬ 
ables, the inclusive spin structure functions gi and g 2 of 
the nucleon are a vital ingredient for this picture (for a re¬ 
view, see my For a complete understanding of the parton 
structure of the nucleon, we need precise and comprehen¬ 
sive data not only for the proton, but also for the neutron. 
Since the two nucleons are isospin partners, one can infer 
(assuming approximate isospin symmetry) the relative con¬ 
tribution from up and down valence quarks as a function of 
momentum fraction x from measurements on protons and 
neutrons. Furthermore, fundamental sum rules concerning 
the difference between proton and neutron structure func¬ 
tions at all values of squared four-momentum transfer 
can be tested experimentally. The isoscalar sum of proton 
and neutron spin structure functions in the Deep Inelastic 
Scattering (DIS) region is particularly sensitive, via pertur¬ 
bative QCD evolution equations IlHl], to the gluon helicity 
distribution inside a longitudinally polarized nucleon. Mo¬ 
ments of structure functions from proton and neutron access 
different matrix elements of local operators within the Op¬ 
erator Product Expansion approach [SHTj. Finally, a better 
understanding of the phenomenon of quark-hadron dual¬ 
ity El [g requires detailed studies of polarized as well as 
unpolarized structure functions of both nucleons in the res¬ 
onance and DIS regions. While suitable free neutron targets 
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do not exist, one can extract spin structure functions for a 
bound neutron using polarized nuclei like and ^He, using 
some prescription to account for Fermi-motion and the ef¬ 
fective polarization of nucleons in nuclei. The results will be 
further affected to some extent by Final State Interaction 
(FSI) effects that are presently unknown. They have been 
estimated to be small in the DIS region [10] but may be 
larger in some part of the kinematic region covered by the 
data reported here. In the following, we quote results for 
the bound neutron without correcting for such FSI effects. 

The CLAS (CEBAF Large Acceptance Spectrometer) col¬ 
laboration at Jefferson Lab has collected a comprehensive 
set of spin structure function data on the proton as well as 
the deuteron over a wide range in « 0.05 — 5 GeV^, 
and over a wide range of final state masses W = 1 — 3 
GeV. A comparable data set has been collected for the neu¬ 
tron, using polarized ^He as an effective neutron target and 
the spectrometers in Jefferson Lab's Ha II A [HHls]. How¬ 
ever, nuclear binding effects have to be accounted for in 
a model-dependent way in order to extract neutron struc¬ 
ture functions from nuclear data. In particular, in the res¬ 
onance region where cross sections and asymmetries may 
vary rapidly with W, Fermi smearing makes the extraction 
of neutron results challenging and somewhat ambiguous. 
For those reasons, neutron data extracted using an inde¬ 
pendent method and a different target, namely deuterium, 
are highly desirable, both to check systematic uncertain¬ 
ties and to more directly access the isoscalar combination 
9i + ff? its moments. Some deuteron data in the res¬ 
onance region exist from the RSS experiment [14], albeit 
over a relatively narrow range in Q^. Many other exper¬ 
iments [15H19] have measured spin structure functions of 
the deuteron in the deep inelastic (DIS) region, W > 2 
GeV and > 1 GeV^, or at small x [20]. Very recently, 
the CLAS collaboration has published precise results from 
the EGl-DVCS run on the proton and the deuteron at the 
highest accessible with Jefferson Lab so far |21j . 

With the experiment presented here (dubbed “EGlb") 
we collected a comprehensive data set on deuteron (^^NDa) 
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targets with nearly equal statistical precision and kinematic 
coverage as on polarized protons (^^NHa). The proton re¬ 
sults will be published separately [22]. In this paper, we 
present our final results for the asymmetry Ai(W,Q^) and 
the spin structure function gi{x,Q^) and its moments for 
the deuteron. The data were obtained in Jefferson Lab's 
Hall B during the time period 2000 - 2001. Previously, a 
much smaller data set on the deuteron was collected with 
CLAS in 1998 [23]. The present data set was taken with 
beam energies of 1.6, 2.5, 4.2 and 5.7 GeV. Preliminary re¬ 
sults from the highest and lowest beam energies have been 
published |24ft26| . The present paper includes, for the first 
time, the full data set collected with CLAS in 2000-2001 on 
the deuteron, including experimental and analysis details. 
We also provide, for the first time, our results for the corre¬ 
sponding (bound) neutron structure functions, based on a 
somewhat model-dependent deconvolution procedure which 
accounts for Fermi motion in the deuteron (27| . 

Our analysis of the deuteron data follows closely that for 
the proton data taken at the same time. Insofar as both 
analyses share the same ingredients and methods, only a 
brief summary is given here - the details will be provided in 
the companion proton paper |22| . However, where the two 
analyses differ, we give all details specific to the deuteron 
in what follows. After a brief summary of formalism and 
theoretical background (Section II), we describe the experi¬ 
mental setup (Section III) and the analysis procedures (Sec¬ 
tion IV). We present the results for all measured and derived 
quantities, as well as models and comparison to theory, in 
Section V, and offer our conclusions in Section VI. 


II. THEORETICAL BACKGROUND 
A. Formalism 


We define the usual kinematic quantities in inclusive 
lepton scattering: Incident (E) and scattered (£’') lep¬ 
ton energy in the lab, scattering angle 0, energy transfer 
V = E — E' and squared four-momentum transfer 

= -q^ =q^ = AEE' sin^ (1) 

The invariant final state mass is 

W = y/M2 + 2Miy-Q2, (2) 


and the Bjorken scaling variable 


^ 2Mv 

in which M is the nucleon mass, 
are also used: 


( 3 ) 


The following variables 


and the virtual photon polarization ratio 

£ = ^1 -l- 2[1 -(- t] tan —^ . (5) 


B. Cross sections and asymmetries 


The observable measured in EGlb is the double spin 
asymmetry 


A^^{iy,Q\E) 


dcrl-'f + dcrlf' 


( 6 ) 


for inclusive electron deuteron scattering with beam and 
target spin parallel (tfi-) or antiparallel (flj.) along the 
beam direction. It depends on the four structure functions 
Ef,Ei,g‘landg^^. Introducing the ratio R of the longitu¬ 
dinal to transverse virtual photon absorption cross sections. 


R=^ 


(J 'J' 


F2 

2xFi 




( 7 ) 


and the variables 


D = 


1 - E'e/E 
lEeR 


and g 


E-E'e' 


we can express A|| as: 



(l + 77)|r + [7(7-7)]|r- 


(8) 

( 9 ) 


Alternatively, the double spin asymmetry A^^ can also be 
interpreted in terms of the two virtual photon asymmetries 
Ai and A 2 : 


A\\ = D[Ai{iy,Q‘^) +gA2{v,Q^)]. 


( 10 ) 


Because of the relative size of the kinematic factors in 
Eqs. [ 9 HT 0 I our data are mostly sensitive to gi or Ai, w hich 
are the main quantities of interest (see Sections |ll C| and 
Given a model or other information for Fi, R and A 2 , 


Ai can be extracted directly from Eq. 10 and gi from 


ffi = + (7-7)^2^ 7^1- (11) 

Our deuteron data are not sensitive enough to A 2 or 32 to 
constrain these quantities; instea d a m odel based on other 
existing data is used (see Section VD|. 


C. Virtual photon absorption asymmetry 

The asymmetry Ai can be interpreted in terms of tran¬ 
sition amplitudes to specific final states (at IT < 2 GeV, 
i.e. in the resonance region) or in terms of the underly¬ 
ing quark helicity distributions (at larger W and Q^). In 


^ In principle, the tensor structure function bi also enters in the de¬ 
nominator, since any realistic polarized target will have a non-zero 
tensor polarization Pzz- However, in our case this is a sub-percent 
correction since Pzz is expected to be less than 0.1 for our tar¬ 
get m and the tensor asymmetry Azz was measured by HERMES 
to be of order 0.01 — 0.02 |29| . 
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the former case, the measured asymmetry Ai at a given 
value of W gives information on the helicity structure of 
the combined resonant and non-resonant contributions to 
the inclusive cross section, which can help to constrain the 
spin-isospin structure of nucleon resonances. 

In the DIS region, Ai{x) can yield information on the 
polarization of the valence quarks at sufficiently large x 
(x > 0.5), where they dominate. In the naive parton model, 
without taking nuclear effects into account, the limit of 
Aid{x) at large x is given as 


to be some duality between structure functions in the res¬ 
onance region (averaged over a suitable range in W) and 
their extrapolated DIS values at the same quark momen¬ 
tum fraction a; or ^ . This correspondence should 

be tested for both nucleon species and for polarized as well 
as unpolarized structure functions to elucidate the underly¬ 
ing dynamics. EGlb data have uniquely suitable kinematic 
coverage stretching from the resonance to the DIS region to 
test whether duality holds for gi. (An initial study of duality 
based on part of the EGlb data has been published |25|.f 


Aid 


Auv + Ady 

Uy dy 


Auy jUy i^dy jUy ) Ady jdy 

1 -h dy/Uy 


( 12 ) 


E. Quasi-elastic scattering 


where Uy, dy are the unpolarized up and down valence quark 
distributions and Auy,Ady are the corresponding helic¬ 
ity distributions. In a S'f7(6)-symmetric, non-relativistic 
quark model |30J, Au/u — 2/3 and Ad/d = —1/3, and 
d/u — 1/2, yielding Aid — 1/3. On the other hand, 
more advanced quark models predict that Aid(x) —?■ 1 as 
a; 1 due to SU( 6 ) symmetry breaking [31]. However, even 
relativistic constituent quark models |32| predict a much 
slower rise towards Ai — 1 than perturbative QGD calcula¬ 
tions |^[M] incorporating helicity conservation. Recently, 
modifications of the pQGD picture to include orbital angular 
momentum [35] have yielded an intermediate approach to¬ 
wards X = 1. Precise measurements of Ai at large x and in 
the DIS region are therefore required for protons, deuterons 
and neutrons to establish the validity of these predictions. 


D. The spin structure function gi 

The structure function gi{x,Q^) contains important in¬ 
formation on the internal spin structure of the nucleon. In 
the DIS limit (large and ly), it encodes the polarized 
Parton Distribution Functions (PDFs) Aq{x) = qt {x) — 
ql {x) for quarks with helicity aligned vs. antialigned with 
the overall (longitudinal) nucleon spin. Its logarithmic 
dependence contains, via the QGD evolution equations Hi¬ 
ll], information on the analogous helicity-dependent gluon 
PDFs AG{x) as well. The deuteron, as an approximate 
isoscalar nucleon target, is particularly sensitive to AG{x), 
given a sufficiently large range in Q^. Jefferson Lab data, 
like those presented in this paper, can serve as a valuable 
anchor point at the lowest possible for NLO fits to ex¬ 
tract Aq{x) and AG{x). 

In the region of lower Q^, additional scaling violations oc¬ 
cur due to higher-twist contributions, leading to correction 
terms proportional to powers of l/Q^. These corrections 
can be extracted from our data since they cover seamlessly 
the transition from ^ 1 GeV^ to the scaling region 
> 1 GeV^. In the kinematic region where iv is also small 
and therefore W < 2 GeV, the structure of 51 is dominated 
by the contributions from nucleon resonances (similarly to 
Ai). 

However, as already observed by Bloom and Gilman [5] for 
the unpolarized proton structure function F 2 , there seems 


The virtual photon asymmetries Ai and A 2 are also de¬ 
fined for elastic scattering off the nucleon and the same 
relationship Eq. applies. One can show that Ai = 1 \n 
this case, and 


A2iQ^) = VR 


GsiQ^) 

V^Gm{Q^)’ 


(13) 


where Ge and Gm are the electric and magnetic Sachs form 
factors of the nucleon. 

One can also extend the definition of gi{x) and g 2 {x) for 
the nucleon to include elastic scattering, x = 1: 

gS‘(^,Q^) = ^^AA-AAs(r-l). (14) 

Z i -F T 

For a bound system like deuterium, one has to consider 
the initial state (Fermi-) motion of the struck nucleons. In 
quasi-elastic inclusive scattering, W < \ GeV, both the 
neutron and the proton contribute (weighed by their elastic 
cross sections). Alternatively, if one detects the struck pro¬ 
ton in addition to the scattered electron with small missing 
four-momentum, the asymmetry i4|| will be close to that on 
a free proton [35]. In both cases, the theoretical asymmetry 
can be calculated with reasonable precision (given a real¬ 
istic deuteron wave function) and therefore the measured 
asymmetry can be used to extract the product of target 
and beam polarization (see below). 


F. Moments 

In addition to the structure function gi{x) itself, its mo¬ 
ments (integrals over x weighted by powers of a:) are of great 
interest. Within the Operator Product Expansion formal¬ 
ism, these moments can be related to local operators [315]. 
They are constrained by several sum rules and can be calcu¬ 
lated directly within lattice QGD or in effective field theories 
like Chiral Perturbation Theory (xPT) pZl|38]- Determin¬ 
ing these moments over a range of allows us to study 
the transition from hadronic degrees of freedom at large dis¬ 
tances (small Q^) to partonic ones at small distances in our 
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description of the nucleon, and to extract higher twist ma¬ 
trix elements that are sensitive to quark-gluon correlations 
in the nucleon. 

The first moment of gi, 

ri(Q^)= [ giix,Q'^)dx, (15) 

Jo 

can be related to the contribution AS of the quark helicities 
to the nucleon spin in the limit of very high . In particular, 
for the average of proton and neutron (the isoscalar nucleon 
approximated by the deuteron) one has 


r?+”(Q2 ^ oo) 
2 


rf = ^{Au + Ad) + ^As. (16) 


understanding of strongly interacting matter at all length 
scales. 

Extending the analysis of low-energy Compton amplitudes 
to higher orders, one can get additional generalized sum 
rules [46]. In particular, one can generalize the forward spin 
polarizability, 70 , to virtual photons: 

70(Q^) = [ x"^ [giix.Q"^)-"f‘^g2{x,Q‘^)]dx. 

W Jo 

( 20 ) 

Once again, this generalized spin polarizability can be cal¬ 
culated using xPT [38] . 


G. From nucleons to the deuteron 


Forming the difference between proton and the neutron 
yields the famous Bjorken sum rule [391140] : 

r?-r” = ^03 = 0.211 (17) 

where = gA = 1.267 ± 0.004 is the neutron axial beta 
decay constant. 

At high but finite Q^, these moments receive logarithmic 
pQCD corrections. At the more modest of our data, 
additional corrections due to higher twist matrix elements 
and proportional to powers of l/Q^ become important: 

/Vf2 

ri(Q^) = A2(Q^)+^^^ [o-2iQ^) + 4d2(Q^) + 4/2((5^)] • • • 

(18) 

Here, /i 2 is the leading twist contribution given by Eq. 
plus pQCD corrections, 02 and d 2 are due to target mass 
corrections and /2 is a twist-4 matrix element that contains 
information on quark-gluon correlations and has been cal¬ 
culated using quark models [41], QCD sum rules (42j and 
other approaches like lattice QCD [43]. 

In addition to the leading first moment, odd-numbered 
higher moments of gi can be defined as T” = 
/p dxx^~^gi{x), n = 3,5,7,.... These moments are dom¬ 
inated by high X (valence quarks) and are thus particularly 
well determined by data in Jefferson Lab kinematics. They 
can also be related to hadronic matrix elements of local op¬ 
erators or evaluated with Lattice QCD methods. The third 
moment Tf is related to the matrix element 02 above. 

In the limit of very small photon virtualities moments 
of spin structure functions can be connected to observ¬ 
ables in Compton scattering. In particular, the first mo¬ 
ment is constrained by the Gerasimov-Drell-Hearn (GDH) 
sum rule [441143 in the limit 0 : 


dTiiQ^) 

dQ^ 


Q2=0 


8M2’ 


(19) 


where n is the anomalous magnetic moment of the nucleon. 
Higher order derivatives at the photon point are, in princi¬ 
ple, calculable via xPT [mi^. Therefore, measuring Ti 
over the whole range in yields a stringent test of our 


Most of the previous discussion is focused on the interpre¬ 
tation of spin structure functions of the nucleon (proton and 
neutron). Where appropriate, we indicate how this interpre¬ 
tation may be modified when the nucleons are embedded in 
deuterium. Here, we want to discuss in more detail how 
the nuclear structure of the deuteron affects the measured 
asymmetries and structure functions. 

In the most simple-minded picture, all observables on the 
deuteron can be considered (cross section weighted) aver¬ 
ages of the corresponding proton and neutron observables. 
Spin observables are further modified by the fact that even 
in a fully polarized deuteron, the nucleon spins are not 100 % 
aligned due to the D-state component of the wave function. 
To first order, this can be corrected by applying a reduction 
factor (1 — 1.5 Pd) to all nucleon spin observables inside 
deuterium [47], with Pd « 4 — 6 % being the D-state prob¬ 
ability (according to the results from recent nucleon-nucleon 
potentials [48]). Taking this factor into account, the spin 
structure functions gf{x) and g^ix) of the deuteron are rea¬ 
sonably well approximated by the average of the proton and 
neutron ones, as long as x is not too large {x < 0.6) and W 
is not in the resonance region (/.e., W > 2 GeV). Moments 
of these structure functions can be considered as relatively 
“safe" since the integration averages over effects like Fermi 
motion Hz). 

In the valence region of moderate to large x and in the res¬ 
onance region, Fermi-smearing due to the intrinsic motion 
of the nucleons inside deuterium as well as nuclear binding 
and FSI become more important, because structure func¬ 
tions vary rapidly in this region with W or x. These binding 
effects can be partially modeled by convoluting the free nu¬ 
cleon structure functions with the momentum distribution 
of nucleons inside deuterium. In our analysis, we use a re¬ 
cent convolution model by Melnitchouk et al. [271 09] that 
properly treats the effects of finite momentum transfer Q^. 

On the other hand, no universal model of the effects of 
FSI over the whole kinematic region covered by our data is 
available; we therefore do not correct for those effects. Sim¬ 
ilarly, potential off-shell effects (due to the negative bind¬ 
ing energy of nucleons inside deuterium), including perhaps 
a modification of the nucleon structure (the EMC effect) 
and non-nucleonic degrees of freedom (mesons [50], AA 
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components [511 f52| and perhaps more exotic quark struc¬ 
tures ^3]) may play a role. Since no universally accepted 
model for these effects exists, we present our results with 
the caveat that they are for bound neutrons only. Given the 
small binding energy (-2.2 MeV) and large average inter¬ 
nucleon distance (of order 4 fm) in deuterium, we expect 
these effects to be significantly smaller than in more tightly 
bound nuclei. However, a comparison with neutron spin 
observables obtained from measurements on ^He can be a 
valuable check on the size of nuclear binding corrections. 
Ultimately, the best approach to extracting free neutron 
information would be to apply the method of spectator tag¬ 
ging (pioneered for unpolarized structure functions in the 
recent “BONuS" experiment p4] at Jefferson Lab). 


III. THE EXPERIMENT 

The EGlb experiment took place at Jefferson Lab over a 
seven month period in 2000-2001. It used the highly polar¬ 
ized (up to 85%) electron beam produced by the Continuous 
Wave Electron Beam Accelerator (CEBAF), with energies 
from 1.6 GeV to nearly 6 GeV and currents of 0.3 to 10 
nA in the experimental Hall B. Detailed descriptions of the 
accelerator and its strained GaAs polarized electron source 
can be found in Refs. |55)J58| . 

The beam polarization was intermittently monitored us¬ 
ing a Mpller polarimeter, and the beam position and inten¬ 
sity distributions were measured with a set of beam mon¬ 
itors. The amount of beam charge delivered to the Hall 
for a given time interval was measured with a Faraday cup 
(FC). The signal from this FC was recorded separately for 
each beam polarization and gated by the data acquisition 
live time. In order to avoid local heating and depolarization, 
the beam was rastered over the face of the target in a spiral 
pattern, using two magnets upstream from the target. 

The target consisted of cells containing samples of po¬ 
larized hydrogen (^^NHa), deuterium (^^NDa), carbon, or 
no solid material (“empty target") that could be alterna¬ 
tively inserted in the beam. These cells were suspended in 
a liquid "^He bath at about 1 K. The target material was 
polarized inside a 5 T solenoidal field along the beam axis, 
using the method of dynamic nuclear polarization (DNP) 
described in [59H61j . The target polarization was monitored 
by an NMR system. Typical values of about 30% deuteron 
polarization along or opposite to the beam direction were 
achieved during the experiment. 

Scattered electrons (and other particles) were detected 
with the CEBAF Large Acceptance Spectrometer (CLAS) 
[62] in Hall B. CLAS employs a toroidal magnetic field 
and several layers of detectors in six identical sectors sur¬ 
rounding the beam axis for an acceptance of nearly 27r in 
azimuth. Electrons were detected in the scattering angle 
range from 8 ° to about 50°. Three regions of drift cham¬ 
bers (DC) [63| determine charged particle trajectories, fol¬ 
lowed by Cherenkov counters (CC) [64] and electromagnetic 
calorimeters (EC) |65| for electron identification, while tim¬ 
ing is provided by a scintillation counter (SC) system [ 66 ] . 


For EGlb, the trigger was optimized for inclusive electrons 
and required a coincidence between signals above threshold 
in the EC and the CC. 

The experimental setup and operation will be described 
in detail in the upcoming companion paper on our proton 
results [ 22 | . 


IV. DATA ANALYSIS 
A. Data set 

Data on the deuteron (ND 3 ) were taken with seven differ¬ 
ent beam energies and two opposite polarities of the CLAS 
torus magnetic field. For positive (-{-) polarity, electrons are 
bent towards the beam line, and for the negative (-) polarity, 
away from it. The in-bending (-f) configuration gives ac¬ 
cess to the largest scattering angles and allows CLAS to run 
with its highest possible luminosity of L = 2 - 10 ^"^ cm^^s“^. 
Therefore, we used this configuration to collect the highest 
points for each beam energy. In the out-bending (-) 
configuration, electrons were detected down to the smallest 
accessible scattering angle of 8 °, extending the data set to 
lower Q^. 

In all, data were collected in 11 specific combina¬ 
tions (1.606-F, 1.606-, 1.723-, 2.561-F. 2.561-, 4.238-F, 
4.238-, 5.615-F, 5.725-F. 5.725-, 5.743-) of beam energy 
(in GeV) and main torus polarity (-F,—), hereby referred 
to as “sets". Sets with similar beam energy comprise four 
groupings with nominal average energies of 1.6, 2.5, 4.2 and 
5.7 GeV. The kinematic coverage of the data for each of the 
4 energy groupings is depicted in Fig. 


B. Data selection 

After following the standard calibration procedures for all 
CLAS detector elements, the raw data were converted into 
a condensed data summary tape (DST) format containing 
track and particle ID information. Quality checks ensured 
that malfunctioning detector components, changes in the 
target and/or potential sources of false asymmetries did not 
contaminate the data. DST files not meeting the minimal 
requirements were eliminated from analysis. 

Event selection criteria were applied to identify scattered 
electrons and to minimize the background from other par¬ 
ticles, primarily tt”. These criteria, based on the signals 
from the CC and the EC, will be discussed in detail in [22] . 
We ascertained that the remaining 7 r“ contamination of our 
electron sample was less than 1 % over the whole kinematic 
range. For this reason, we assign a 1% systematic uncer¬ 
tainty on our extracted asymmetries as an upper limit for 
any remaining pion contamination effect. 

For the determination of the product of beam and tar¬ 
get polarization [PbPt, see below) as well as kinematic cor¬ 
rections, we also required a sample of quasi-elastic {e,e'p) 
events. We selected ep coincidences through a timing cut 
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FIG. 1. (Color Online) Kinematic coverage in vs. x for 
each of the 4 main electron beam energy groupings used in 
the EGlb experiment. The solid and dotted lines denote the 
W = 1.08 and W = 2.00 GeV thresholds, respectively. The 
coverage for proton (NH 3 ) and deuterium (ND 3 ) targets was 
nearly identical. 


xIO'® 



<P^i - (Ppr (degrees) 


FIG. 2. (Color Online) Distribution of quasielastic d{e, e'p) 
events versus the angle (f) between the azimuth of the scat¬ 
tered electron and the azimuth of the observed proton. The 
background due to nitrogen, liquid ^He and various foils is 
strongly suppressed by the cuts described in the text, lead¬ 
ing to a relatively clean signal from the deuteron component 
(solid line) of the target. A final cut is applied from (j) = 177° 
to 183°. 


of ±0.8 ns on the difference between the reconstructed elec¬ 
tron and proton vertex time. Quasi-elastic events were se¬ 
lected through cuts on W, 0.89 GeV < W < 1.01 GeV, 
missing energy (of the unobserved nuclear remnant) of 
< 0.08 GeV (kinetic), and on the difference between the 
polar (|A0| < 2°) and azimuthal (|A(?!)| < 3°) angles of 
the detected proton and the reconstructed direction of the 


virtual photon. These cuts were optimized to include most 
of the ep coincidences from quasi-elastic scattering on the 
deuteron, while the contribution from the other target com¬ 
ponents (nitrogen, "^He and foils) was much suppressed due 
to the wider nucleon momentum distributions in these nuclei 
(see Fig. [^. 

C. Event corrections 

The track information for particles in the DSTs is based 
on an ideal detector and has to be corrected for various 
effects from detector materials and imperfections. Among 
other corrections, energy loss due to ionization in the target 
(both for the incoming and the scattered electron), multiple 
scattering angle deviations (compared to the average vertex 
of all particles in an event), and known deviations of the 
target magnetic field from the ideal version implemented in 
the reconstruction software were used to correct each track 
within an event. 

The reconstruction software also assumes that a track 
originates on the nominal central axis (a; = y = 0) of CLAS. 
In reality, the beam is rastered over a circle of about 1.5 cm 
diameter, whose center is typically offset by a few mm from 
the nominal axis. Since the raster position can be inferred 
from the currents in the raster magnets, the reconstructed 
vertex was corrected for this offset. 

The position and orientation of the drift chambers in 
space and the detailed three-dimensional shape of the 
torus magnetic field are not known with absolute preci¬ 
sion; an empirical parametrization of their deviations from 
the ideal detector was obtained from a fit to data from 
the companion experiment on the proton [22]. We used 
four-momentum conservation in fully exclusive events like 
H(e, e'p) and H(e, e'pTr'^TT~) to optimize the fit parameters. 
This parametrized correction for particle momenta and scat¬ 
tering angles was then applied to each track. The resulting 
improvement of the resolution in the missing mass W is 
shown in Fig. 

A final correction was applied to the integrated beam 
charge measured by the Faraday Gup, to account for beam 
loss between the target and the FG due to multiple scatter¬ 
ing and due to dispersion by the target magnetic field. 

D. From raw to physics asymmetries 

For each combination of beam energy, torus polarity and 
target polarization, electron tracks were sorted by kine¬ 
matic bins and were counted separately for positive (iV+) 
and negative {N~) beam helicity, where “±" refers to a 
beam helicity antiparallel to the direction of the target 
polarization. These counts were normalized to the cor¬ 
responding integrated Faraday charges, = N^/FC^. 
Only events coming from complete pairs of “beam buckets" 
with opposite helicity were counted to avoid false asym¬ 
metries; we also ascertained that, after averaging over all 
target polarizations, the residual beam charge asymmetry 
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FIG. 3. (Color Online) Missing mass W before (red-hollow) 
and after (blue-solid) the kinematic corrections for the 4.238-1- 
data set for NH 3 (top) and ND 3 (bottom) targets. The correc¬ 
tions decreased the distribution width and centered the mean 
value of the (quasi-)elastic peak on the nucleon mass. 

(FC+ - FC-)/(FC+ + FC-) was less than IQ-^. These 
normalized counts were used to form the raw asymmetry 

_ n+-n- 

■^raw — I , 

n+ + n~ 

in each kinematic bin. This raw asymmetry was then con¬ 
verted to the desired physics asymmetry zl|| (Eq.[^ by ap¬ 
plying a series of corrections which we now discuss in se¬ 
quence. 


1. Dilution factor 

The dilution factor Fjjp = Ud/riA is defined as the ratio 
of events from polarizable nuclei of interest (here, deuterons 
bound in ammonia, Ud) to those from all components of the 
full ammonia target (n^). It is calculated directly from the 
radiated cross-sections on all components of the target. In 
terms of densities (p), material thicknesses {i) and cross- 




FIG. 4. (Color Online) Dilution factors as a function of 
W, shown at four different beam energies (1.6-1- (top left), 
2.5— (top right), 4.2— (bottom left) and 5.7— (bottom right)). 
The results from our standard method (using cross section 
models, see text) are shown as blue lines, while the results 
from the data-based method (see text) are shown as the red 
data points. 

sections per nucleon (a), 

nd oc ^pA^AfJd (22) 

and 

nA oc PAl^AlCtAl + PK^KCTk 
6 15 

+ PA^Ai^Ctd + + PHe{L — (■A)t^He, (23) 

with the subscripts A, Al, K, N, and He denoting deuter- 
ated ammonia (^®ND 3 ), aluminum foil, kapton foil, nitro¬ 
gen (^®N) and helium ("^He), respectively. The acceptance- 
dependent proportionality constant is identical in both of 
the above relations for a given kinematic bin. Inclusive 
scattering data from the empty (LHe) and targets were 
analyzed to determine the total target cell length (L) and 
effective ND 3 thickness {£a) using similar equations. 

The required cross-sections were calculated from a fit to 
world data for Fi and F 2 for protons and neutrons, using 
a Fermi-convolution model to fit inclusive scattering data 
on nuclear targets, including EGlb data from solid 
and empty (LHe) targets The nuclear EMC ef¬ 

fect was parametrized using SLAC data [69j . Radiative cor¬ 
rections used the treatment of Mo and Tsai IZQ]; external 
Bremsstrahlung probabilities incorporated all material thick¬ 
nesses in CLAS from the target vertex through the inner 
layer of the DC. 

Dilution factors Fjyp were calculated for each data set 
and used to correct the raw asymmetry, 

( 24 ) 

^DF 
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to get the undiluted asymmetry due to deuterons in the 
target. We checked our results for Fdf from the “stan¬ 
dard method” described above against a previously devel¬ 
oped "data-based method" [26l ITT] that uses a simple 
model of neutron/proton cross-section ratios to express the 
background in the ammonia target in terms of the counts 
from carbon and empty targets. Values of L and varied 
by less than 2% between the two methods. Figure]^ shows 
the result from both methods for four kinematic bins. For 
the inelastic data, W > 1.1 GeV, the dilution factors from 
the cross-section based standard method were more pre¬ 
cise and were used to correct the raw asymmetries. We 
used the data-based method only in the quasi-elastic region 
W < 1.08 GeV (for the determination of beam and tar¬ 
get polarization in one case, see below) and to subtract the 
background from exclusive d{e,e'p)n events (see Fig. |^. 
This is because finite detector resolution effects (which are 
not included in the cross section model) significantly affect 
the shape of sharply peaked spectra in the quasi-elastic re¬ 
gion, making the data-driven method more reliable. 

The densities and thicknesses of all target materials were 
varied within their known tolerances to determine system¬ 
atic uncertainties; only the variations of pcic and pHe had 
any significant (>1%) effect on Fdf- Uncertainties due to 
the cross-section model were estimated by the comparison 
of Ff)f to a third-degree polynomial fit to the data-based 
dilution factors determined by the alternate method. 


2. Beam and target polarizations (PbPt) 

The second major factor to consider when extracting the 
physics asymmetry A|| is the product of beam and target 
polarization by which the measured asymmetry must be di¬ 
vided. 

Because NMR measurements provided accurate target 
polarization measurements only near the edge of the target 
cell [72) (which was not uniformly exposed to the beam), we 
determined the polarization product PbPt directly from our 
data, using quasi-elastic (i(e, e'p)n and (in one case) d(e, e') 
events. Here, we made use of the fact that the theoretical 
asymmetry in this case depends only on the electromagnetic 
form factors of the proton and the neutron, see Section m 
which are well-known [73], giving us reliable predictions of 
A||. After correcting for the (relatively smaller) dilution 
of this asymmetry from non-deuterium components of the 
target, we can directly divide the measured A|| by the the¬ 
oretical one to extract PbPt- 


PbPt 


^QE 

meas 

Fdf 


(25) 



- i - 

-A- 

-J 

.. T"" 

■ ■T' 



_ 


1 
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f 
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n ^ 1. 
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FIG. 5. (Color Online) PbPt values for the 2.5 GeV in¬ 
bending data sets. The plot shows the resulting PbPt values 
extracted independently for each bin with available data. 
The results from the exclusive (blue filled symbols) and the 
inclusive (red open symbols) methods are shown. The corre¬ 
sponding constants fit to the data are also shown as lines: the 
solid blue line is for the exclusive and the dashed red line is 
for the inclusive method. 


yields a smaller statistical uncertainty, it has greater system¬ 
atic uncertainty because of larger background contributions; 
therefore, a systematic uncertainty of 10% was assigned to 
this particular PbPt value. 

For all other configurations, we determined PbPt using 
exclusive d{e,e'p)n events within the cuts listed in Sec¬ 
tion [ly[^ which have very little background from nuclear 
target components (see Fig.[^. We used a detailed Monte 
Carlo simulation, including Fermi motion of the proton in¬ 
side the deuteron, to calculate the theoretical asymmetry. 
For both methods, the nuclear background was determined 
using the data-driven method mentioned in Section p^V D 1| 
As a cross check, we compared these results to the values 
derived from inclusive quasi-elastic scattering, and found 
them generally to be consistent within the statistical uncer¬ 
tainty. 

The derived PbPt values were checked for consistency 
across for each beam energy, torus current and target 
polarization direction. Sample PbPt values across for 
2 beam energies are shown in Figure]^ Across all beam 
energies, PbPt values ranged from 0.1 to 0.28, with most 
values between 0.15 and 0.25. We varied each of the values 
of PbPt individually by the larger of one (statistical) stan¬ 
dard deviation and the difference between the exclusive and 
inclusive results to assess the systematic uncertainty of all 
physics quantities due to PbPt- 


3. Polarized nitrogen and target contamination corrections 


We used the value for PbPt obtained from inclusive quasi¬ 
elastic events only in one case, for the 1.6 -1.7 GeV outbend- 
ing configuration runs. In that case, too few of the protons 
from d{e, e'p)n were detected in GLAS for a reliable deter¬ 
mination of PbPt- We used a cut of 0.89 GeV < W < 
1.01 GeV to define quasi-elastic events. While this method 


Apart from the dilution of the measured asymmetry by 
nucleons embedded in nitrogen, helium and other target 
materials (Section IVD1|, there are additional small modi¬ 
fications of this asymmetry due to polarized target nucleons 
outside of deuterium. 

First, it is well-known that the nuclei in the ammo- 
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nia molecules become somewhat polarized as well. Equal 
Spin Temperature (EST) theory predicts the polarization 
ratio between two spin-interacting nuclear species in a ho¬ 
mogenous medium as the ratio of their magnetic moments: 

« 1115]S!/ii 2 u. However, experimentally, it was 
found that the polarization is somewhat smaller than 
that [74| . Using a simple shell model description [75] of the 
nucleus, this polarization is carried by a single proton 
in the lpi /2 shell, which means that this proton is spin- 
polarized to -33% of the nucleus. The measured magnetic 
moment of suggests a somewhat smaller spin polar¬ 
ization, so that the overall contribution from nitrogen to 
the measured asymmetry can be approximated by that of 
a bound proton with polarization p^ound between 8% and 
16% of the deuteron polarization. Accordingly, we sub¬ 
tracted a correction of 1/3 x x Apcrp°“”‘^/crd « 

(0.026 ± 0.014)Ap from the measured asymmetry, where 
the factor 1/3 accounts for the three deuteron nuclei per 
nitrogen nucleus in ammonia. 

A second contamination to the measured asymmetry 
comes from isotopic impurities of the deuterated ammonia, 
with some deuterons replaced by protons. Typical contam¬ 
inations quoted in the literature [15] are around 1.5%. We 
did a careful study [T^ that showed a contamination of 
up to about 3.5% during EGl (which was included in the 
dilution factor); however, according to this study at most 
one-half of these extra protons were polarized (the remain¬ 
der are presumably bound in molecules like H 2 O and are 
unpolarized). The degree of polarization of these protons 
can be estimated as Pp/Pd ~ 1-2 — 1.5, again according to 
EST and empirical evidence [75] . The net effect is an ad¬ 
ditional term proportional to Ap that has to be subtracted 
from the measured asymmetry. The total correction for 
bound and free polarized protons in the target is between 
0.027Ap and O.OSlAp. We took the median of this range 
to correct our data (using a model of the asymmetry Ap 
based on our proton results |^) and 1/2 of its spread to 
estimate systematic uncertainties. An additional correction 
due to the very small contribution of nuclei (less than 
2 % of our ammonia sample) was too small to be applied 
but was included in the overall systematic uncertainty. 

Quasi-elastic d{e,e'p)n events are also affected by the 
various target contaminations discussed above. We applied 
a corresponding correction to our extraction of PbPt (Sec¬ 
tion 


IVD2I. 


4. Other background subtractions 

Dalitz decay of neutral pions m and Bethe-Heitler pro¬ 
cesses [78] can produce e+e” pairs at or near the vertex, 
contaminating the inclusive e~ spectrum. This contamina¬ 
tion was at most a few percent of the data rate (at high 
W) and was measured by comparing positron and electron 
rates for runs with opposite torus polarity. We also mea¬ 
sured the positron asymmetry and found it consistent with 
zero. We subtracted this pair-symmetric background using 
the measured rate and assuming zero asymmetry. To esti¬ 


mate the corresponding systematic uncertainty, we instead 
applied a correction assuming a constant positron asymme¬ 
try within the range of values we measured. We also used 
the change in the correction after varying the rate within its 
uncertainty as a second contribution to the overall system¬ 
atic uncertainty for this background. 



1.5 2 2.5 

W (GeV) 


FIG. 6. (Color Online) Representative results for the fully 
corrected double-spin asymmetry A|| versus final state in¬ 
variant mass W for three different bins and beam en¬ 
ergies. The red-solid li ne re presents our model parametriza- 
tion of A|| (see Section VD l. The dashed blue lines represent 
the model including radiative effects. The difference between 
those lines corresponds to the magnitude of radiative correc¬ 
tions applied. The error bars reflect statistical uncertainties 
while the shaded bands at the bottom of each plot represent 
the total systematic uncertainties. 


5. Radiative corrections 

Radiative corrections to the measured asymmetries A|| 
were computed using the program RCSLACPOL, which was 
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developed at SLAC for the spin structure function experi¬ 
ment E143 ^9j. Polarization-dependent internal and exter¬ 
nal corrections were calculated according to the prescrip¬ 
tions in Ref. izg and Ref. [^, respectively. 

We compared the calculated double spin asymmetry with 
radiative effects turned on, Ar, to the Born asymmetry, Ab, 
calculated with the same models (see Section VD). We 
determined parameters fuc and Arc for each kinematic 
bin, allowing us to write the Born asymmetry as 


Most sources of systematic uncertainties have been dis¬ 
cussed above. These sources include kinematic shifts, bin 
averaging, target parameters (radiative corrections), nuclear 
dilution model, structure function models, PbPt uncertainty 
for each individual data set, and background contamina¬ 
tions. The relative magnitudes of these various contribu¬ 
tions to the systematic uncertainty, for the case of the ratio 
gijFi, are listed in Table [I] The results shown in the next 
section incorporate these systematic uncertainties. 


Ab = -jA- + Arc^ 

IRC 


(26) V. RESULTS AND COMPARISON TO THEORY 


where fuc is a radiative dilution factor accounting for the 
count rate fraction from the elastic and quasi-elastic tail 
within a given bin. This correction was then applied to all 
data. Figure shows a few examples for the magnitude 
of the correction, together with the final data for the Born 
asymmetry A||. 

Systematic uncertainties on these corrections were es¬ 
timated by running RCSLACPOL for a range of reason¬ 
able variations of the models for F 2 , R, Ai and A 2 (see 
Section VD) and for different target thicknesses and cell 
lengths, £a and L. The changes due to each variation were 
added in quadrature and the square root of the sum was 
taken as the systematic uncertainty on radiative effects. 


6. Systematic uncertainties 

Estimation of systematic uncertainties on each of the 
observables discussed in the following section was done 
by varying a particular input parameter, model or analysis 
method, rerunning the analysis, and recording the differ¬ 
ence in output for each of the final asymmetries, structure 
functions and their moments. Final systematic uncertain¬ 
ties attributable to each altered quantity were then added 
in quadrature to estimate the total uncertainty. Note that 
for each quantity of interest (Ai, gi, Ti) the systematic un¬ 
certainty was calculated by this same method (instead of 
propagating it from other quantities), therefore ensuring 
that all correlations in these sources were properly taken 
into account. 


Systematic uncertainty 

Typical range (in % of gi/Fi) 

Pion and e^e“ contamination 

0.0% - 1.0% 

Dilution Factor 

1.8% - 2.7% 

Radiative corrections 

3.5% - 5.7% 

PbPt uncertainty 

6% - 22% 

Model uncertainties 

2.0% - 5.0% 

Polarized Background 

1.0% - 1.7% 

Total 

10% - 23% 


TABLE I. Table of typical magnitudes for various systematic 
uncertainties. 


A. Results for Ai + rjA 2 



1.2 1.4 1.6 1.8 2 

W (GeV) 


FIG. 7. (Color Online) Representative values for the double¬ 
spin asymmetry Ai + rjA 2 versus final state invariant mass 
W . The top panel is for 0.16 GeV^ < <3^ < 0.19 GeV^ (1.6 
GeV data) and the bottom panel for 0.45 GeV^ < <3^ < 0.54 
GeV^ (2.5 GeV data). The red-solid line re prese nts our model 
parametrization of Ai -|- r]A 2 (see Section VDi. The shaded 
band at the bottom (green) is the total systematic uncer¬ 
tainty. The individual contributions are offset vertically, from 
top to bottom: pion and pair symmetric contamination (-0.4; 
barely visible); dilution factor (-0.6); PbPt (-0.8); models plus 
radiative corrections (-1.0); and polarized background (-1.2). 


In this section, we present our final results for all quan¬ 
tities of interest: Ai, gi and moments for the deuteron 
and the bound neutron. As a first step, we divide the fully 
corrected Born asymmetry A|| by the depolarization factor 
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D (Eq. 1^ to extract the combination Ai + 77^2 for each 
bin in W and and each beam energy. Results for sim¬ 
ilar beam energies (e.g., 1.6 and 1.7 GeV) and inbending 
and outbending torus polarization are combined into aver¬ 
aged values for four nominal energies (1.6 GeV, 2.5 GeV, 
4.2 GeV and 5.7 GeV), weighted by their statistical preci¬ 
sion. We checked that in each case, the data sets that we 
combined agree with each other within statistical and sys¬ 
tematic uncertainties. Figures |^and[^show the results for 
Ai -F 77 A 2 for selected bins and for each of the four stan¬ 
dard energies. The systematic uncertainties from different 
contributing sources are also shown as shaded bands at the 
bottom of each plot. For most kinematics, the largest con¬ 
tribution to the systematic uncertainty is due to the beam 
and target polarization, with some contribution from the 
dilution factor and radiative corrections. We note that our 
data for all 4 beam energies are well described by our model 
(see Section VD| as indicated by the red solid line. 

Our results for Ai -F 77 A 2 have the least theoretical bias 
from unmeasured structure functions like A 2 and Fi, and 
are therefore the preferred choice for NLO fits that will in¬ 
clude our data in the high-Q^, high-lV region, like the fit 
by the JAM collaboration [BO]. They can be found in the 
GLAS database [81] and in the Supplemental Material |82l 
for this paper. 



kV (GeV) 

FIG. 8. Same as Fig.[7]except for two higher beam energies. 
The top panel is for 0.64 GeV^ < <3^ < 0.77 GeV^ (4.2 GeV 
data) and the bottom panel for 1.1 GeV^ < < 1-3 GeV^ 

(5.7 GeV data). 


B. The virtual photon Asymmetry Ai 



1.5 2 2.5 

W (GeV) 


FIG. 9. (Color Online) Virtnal photon asymmetry Ai for the 
deuteron versus W for a few bins: 0.16 GeV^ < < 0.19 

GeV^ (top), 0.45 GeV^ < < 0.54 GeV^ (middle) and 1.1 

GeV^ < < 1.3 GeV^ (bottom). The statistical uncertain¬ 

ties are indicated by error bars, while the total systematic un¬ 
certainties are indicated by the shaded band at the bottom. 
Again, the individual contributions are shown separately as 
offset bands: pion and pair symmetric contamination (-0.4); 
dilntion factor (-0.6); Pi,Pt (-0.8); models pins radiative cor¬ 
rections (-1.0); and polarized backgronnd (-1.2). 


Once Ai -F 77 A 2 is calculated, we can extract the virtual 
phot on a symmetry Ai, by using a model for A 2 (see Sec¬ 
tion VD). Since Ai depends only on W and Q^, we can 


combine the results from all beam energies at this stage, 
again weighted by statistical uncertainties. Figureshows 
Ai(W) for three representative bins together with dif- 
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ferent sources of systematic uncertainties. The uncertainty 
on A 2 (included in the band at -1.0) is the dominant contri¬ 
bution to the overall systematic uncertainty (shaded band 
at the bottom of each panel). 

Figures 10 and 11 show Ai versus W for all bins 
in our kinematic coverage, as well as existing data from 
SLAC E143 [ 15 ] |83] and from the Jefferson Lab RSS ex¬ 
periment [m m\- Gaps are due to a lack of kinematic 
coverage between the different beam energies. Data points 
with very large statistical or systematic uncertainties were 
omitted from these plots. 


negative or close to zero, which proves that the A 3/2 transi¬ 
tion amplitude is dominant in this region as expected from 
exclusive pion production. As we go to higher values of 
and W, the transition amplitude Aij^ leading to resonances 
such as A^(1520)3/2” and iV(1535)l/2“ becomes domi¬ 
nant, as expected from pQCD. ktW>2 GeV and larger 
, the asymmetry continues smoothly from the resonance 
region into the DIS region where it has been measured by 
previous experiments to be positive, due to the larger con¬ 
tribution from the proton (with > 0 throughout the 
measured x range in the DIS region). 


K 



K 



W (GeV) 


FIG. 10. (Color Online) Ai for the deuteron versus W for our 
14 lowest bins. Total systematic uncertainties are shown 
as shaded area at the bottom of each plot. Our parametrized 
model is also shown as a red line on each plot. Only the data 
points with tJatat < 0.3 and (Taya < 0.2 are plotted. In addition, 
we also show data from SLAC E143 | 15[ I83 | (open-magenta 
squares). 


At all but the highest Q^, the effect of the A(1232)3/2+ 
resonance is clearly visible in the strongly negative values 
of Ai, due to the dominance of the A 3/2 transition to this 
resonance. At our lowest Q^, the asymmetry is in general 


FIG. 11. (Color Online) Continuation of Fig. 10 for the 
remaining bins. In addition to our data and the SLAC 
data (see above), we also show the data from the Jefferson 
Lab RSS experiment [141184] (blue open circles). 


This trend becomes more apparent if we integrate our 
data on Ai over the full measured DIS range with W > 2 
GeV and > 1 GeV^ and plot it as a function of the scal¬ 
ing variable x. The behavior of Ai{x) at large x is of high 
interest to test various models inspired by QCD, as outlined 
in Section Figure [T^ shows this quantity from EGlb 
together with world data and various models. We note that 
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FIG. 12. (Color Online) Af versus x in the DIS region 
> 1 GeV^ and IF > 2 GeV) from EGlb and several 
other experiments: EGl-dvcs at Jefferson Lab [21], SMC at 
CERN [HI], E143 and E155 at SLAG [Il[l2llH3] and HER¬ 
MES at DESY [TS]. Statistical uncertainties are indicated by 
error bars, and EGlb systematic uncertainties by the shaded 
band at the bottom. Various theoretical predictions and 
parametrizations are shown as lines and shaded band, and 
are discussed in the text. 


our data lie somewhat below most of the world data, which 
is partially explained by the fact that at each point in x, they 
have the lowest average of all the experiments shown, 
implying a more significant impact of scaling violations due 
to higher twist effects. In particular, the new results from 
EGl-dvcs [2T| (also shown in Fig. 12) are for 5.9 GeV beam 
energy and scattering angles above 18°, while our data av¬ 
erage over 5.7 and 4.2 GeV and scattering angles down to 
8°. In addition, systematic differences exist between these 
two most precise data sets due to the target polarization, di¬ 
lution factor, and the different impact from required model 
input for R and A 2 at different kinematics. The correspond¬ 
ing systematic uncertainties are indicated for EGlb by the 
shaded band at the bottom of the plot. 

We also show various predictions based on expectations 
about the asymptotic value for Af in the limit x —?► 1 (see 
Section II C). The prediction from a SU(6)-symmetric quark 
model is a constant value of 1/3 for Af and is indicated 
by the short horizontal line at the right-hand edge of the 
plot. A more advanced quark model including hyperfine 
perturbation through one-gluon exchange [32] yields a range 
of possible behaviors at high x, as indicated by the shaded 
(light blue) band. Two different curves (labeled BBS) are 
based on pQCD models; one under the assumption of pure 
quark-hadron helicity conservation m and a second one 
including the effect of a possible non-zero orbital angular 
momentum (BBS-fOAM [35]). Finally, we show two recent 
NLO parametrizations of the world data (by Soffer [86| and 


by Leader, Stamenov and Sidorov - LSS [87]). 

We note that, on average, the world data including our 
own indicate a rise of Af beyond the SU(6) limit at very 
large x, but much slower than expected from pQCD without 
the inclusion of orbital angular momenta. Taking a possible 
dependence and systematic uncertainties into account, 
our data agree best with the model including orbital angular 
momenta [35] and are also compatible with the lower edge 
of the range of predictions from the hyperfine-perturbed 
quark model [32]. Overall, no firm conclusion can be drawn 
yet about the transition of the down quark polarization from 
negative values below x « 0.5 to the limit of -|-1 expected 
from pQCD. A similar conclusion comes from measurements 
on 3He [12] [88]. 


C. The spin structure function gi 


In addition to extracting Ai, we can also use the mea¬ 
sured asymmetry A\\ to extract the spin structure function 
gf according to Eq.l^ As a first step, we extract the ratio 
gf/Ff which is less sensitive to various model inputs. Fig¬ 
ure]^ shows the resulting data, plotted for several x bins 
(all with a bin width of Ax = 0.05) versus the photon vir¬ 
tuality Q^. Again, we also show world data for the same 
quantity. Our data agree reasonably well with those from 
E143 [151183| within statistical uncertainties, but are some¬ 
what lower than the very precise data from the recently 
published follow-on experiment EGl-dvcs [21]. However, 
the difference is consistent with the (largely uncorrelated) 
systematic uncertainties of both experiments. The de¬ 
pendence at lower reflects the effect of nucleon reso¬ 
nances at IF < 2 GeV, while beyond this limit (indicated 
by arrows on the x-axis) this dependence is mild but still ris¬ 
ing, indicating a smooth but not necessarily fast transition 
to the scaling region. We indicate the results for gifFi at 
= 5 GeV^ from a recent NLO fit of the world data ^7| 
for comparison. 

We then use models for the unpolarized structure function 
Fi (see next section) to convert these ratios to gi. The 
results for the product xgf versu s Bj orken x for each of 
our bins are presented in Fig. 14 together with world 
data. The red curve on each plot comes from our model. 
At low , gi is strongly affected by resonance structures, 
in particular the A(1232) again being the most prominent 
one, making gi negative in this region. When we go to 
higher Q^, the effect of the resonances diminishes and gi 
approaches the smooth DIS curve also shown in Fig. 0 as 
blue dashed line. This can be interpreted as a sign that 
quark-hadron duality begins to work at these larger > 
1.0 GeV^. However, in the A(1232) region, the data fall 
noticeably below the blue line even at as high as « 1 
GeVT 

In the DIS region {W > 2 GeV and > 1 GeV^), gf{x) 
can be used to extract information on the quark helicity con¬ 
tributions to the nucleon spin (see Section pTPj ). Comparing 
our data to the higher data from COMPASS [2^ one 
can extract information on the gluon polarization through 
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FIG. 13. (Color Online) The ratio gi/F\ for the deuteron 
versus and for various bins in the Bjorken variable x, to¬ 
gether with our model shown as the red line on each plot. All 
data are corrected by our model to center them on the middle 
of each x bin. The shaded area at the bottom of each plot 
represents the systematic uncertainty. Published world data 
are shown as open-magenta squares (E143 |15[ 183) 1 and open 
blue triangles (EGl-dvcs [^). Arrows on the x-axis indicate 
the limit W = 2 GeV. The horizontal arrows on the r.h.s of 
the right panel indicate the results for g\ /F\ of a recent NLO 
analysis of world data |87] for our bin centers and = 5 
GeV^. 


DGLAP evolution. Including our data for somewhat lower 
Q^, higher twist modifications of the polarized PDFs can 
be constrained. Our data are available for such PDF fits, 
similar to recent fits by the JAM collaboration [80] and by 
Leader et al. m. as well as for future tests of duality. 


gi (see next section). For the deuteron case in particular, 
we need models for both the proton and the neutron, as 
well as a prescription for Fermi-smearing. 

We will describe our fit in detail in Ref. [22]- Our ap¬ 
proach to Fermi-smearing is explained in Section [IT^ Here, 
we just summarize our sources of data for the fits to A 2 
and Al for the proton and the neutron. For the un¬ 
polarized structure functions Ff’^ and we used a 

recent parametrization of the world data by Bosted and 
Christy [671 189] . This parametrization fits both DIS and 
resonance-region data with an average precision of 2-5%, 
including Jefferson Lab Hall C data on the proton and the 
deuteron with very similar kinematics to ours. Systematic 
uncertainties due to these models were calculated by vary¬ 
ing either Fi or R by the average uncertainty of the fit and 
recalculating all quantities of interest. 

For the asymmetries in the region W > 2 GeV, we de¬ 
veloped our own phenomenological fit to the world data, 
including all DIS results from SLAG, HERA, CERN and Jef¬ 
ferson Lab (see Ref. [1] for a complete list). In the reso¬ 
nance region, we added data from EGla |^|7T] in Hall B, 
RSS [14] in Hall C and MIT-Bates [90]. We also used the 
data reported here and in [22] and iterated the fit after re¬ 
extracting our data using the updated models. The proton 
asymmetries were fit first, followed by a fit to the neutron Ai 
and A 2 . For this second part, we used the rich data set col¬ 
lected on ^He at Jefferson Lab (Hall A) |1 II[T2l [8811911 [92] . 
SLAG |93H96[ . and HERMES [IS] [^, as well as the world 
data on the deuteron, including our own. The goodness 
of the fit (x^) was calculated by comparing the fit func¬ 
tions for neutron asymmetries directly with neutron results 
extracted from ^He data, as well as comparing the convolu¬ 
tion of our proton and neutron models with corresponding 
deuteron data. To anchor our fit of Ai at the photon point, 
we used data from ELSA and MAMI (see, e.g., the summary 
by Helbing [98]). As a result, we achieved a consistent fit 
of proton, deuteron and neutron data over a wide kinematic 
range, far exceeding our own kinematic coverage. The over¬ 
all x^ for the fit was 2451 for 3225 degrees of freedom. 

Our fit results are shown as curves on most of the plots in 
this section, and they are generally in very good agreement 
with the existing data. We developed alternative model fits 
representing the uncertainty of our fit results in all cases 
and estimated the systematic uncertainties on all extracted 
quantities due to model uncertainties by replacing the stan¬ 
dard fits, one by one, with these alternatives. 


D. Models 


E. Moments of gi 


To extract the physics quantities discussed above from 
our data on A||, we require models both for the unpolarized 
structure functions Fi and F 2 (or, equivalently, Fi and R), 
as well as for the asymmetry A 2 . These models (plus a 
model for the asymmetry Ai ) are a lso needed to evaluate 
radiative corrections (Section IV D 51 and to extrapolate our 
data to small x, for the purpose of evaluating moments of 


From our data, we determined several moments of spin 
structure functions. We evaluated those moments for each 
of our standard bins in two parts. For W regions where 
we have good data (with reasonably small statistical uncer¬ 
tainties), we summed directly over these data (binned in 10 
MeV bins in W), multiplied by the corresponding bin width 
in X and the required power of x. We avoided the region 
below W = 1.15 GeV, where radiative effects and the quasi- 
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FIG. 14. (Color Online) The product xgi versus x for all bins, together with our model (red lines). The shaded area at the 
bottom of each plot represents the systematic uncertainty. The corresponding DIS parametrization for — 10 GeV^ is also 
shown (blue dashed lines). World data are shown for Hermes [18] (red circles), SLAG E143 [151183j (open-magenta squares), 
SLAG E155 [T^ (magenta inverted triangles), RSS |14[I84 | (blue circles), and EGl-dvcs (cyan triangles). 


elastic contribution overwhelm the data. The upper end of 
the integration range can go up to VF = 3 GeV, depending 
on the bin. The resulting values of the integral over 
the kinematic region covered by our data are shown as the 
open (blue) circles in Fig. 15 and the properly propagated 
systematic uncertainty in the measured region is shown as 
the light beige band. Note that all moments are calculated 
per nucleon (i.e., divided by 2 for the two nucleons in deu¬ 


terium), following common practice. However, we do not 
correct for the deuteron D-state or any other nuclear effects. 

We integrate our model for gf (without any quasi-elastic 
contributions) over the region 1.08 GeV < W < 1.15 GeV 
in order to estimate this small part of the full moment]^ 


^ We exclude the (quasi-)elastic region W < 1.08 GeV, following 
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FIG. 15. (Color Online) Fi for the deuteron versus from 
our data only (hollow blue circles) and from data plus model 
(full blue circles), including the extrapolation to the unmea¬ 
sured kinematics. The left-hand side shows the full range 
(leaving out our data for < 0.3 GeV^, to avoid clutter) 
and the right-hand side focuses on the small-Q^ region. The 
systematic uncertainty is shown at the bottom of the plot, 
for data only (light beige shaded area in the foreground) 
and for combined data and model (blue shade in the back¬ 
ground). Corresponding results from SLAG E143 |T^ [83|. 
HERMES [18] and EGl-dvcs [21] are shown, as well as sev¬ 
eral predictions (explained in the text). 


Occasionally, there are gaps in our W coverage from differ¬ 
ent beam energies, especially at low (see, e.g.. Fig. 101. 
These gaps are also filled by integrating the model instead. 
Finally, we integrate the model from the lower x limit of 
our highest W bin (for each Q^) down to a: = 0.001. This 
contribution becomes most important at high and for 
the lowest (first) moment. We limit ourselves to this mini¬ 
mum X value because there are no reliable data at lower x, 
and our model becomes unconstrained and rather uncertain 
below X = 0.001. While it is likely that there is no signifi¬ 
cant contribution below this limit we prefer to quote our 
results as moments from x = 0.001 to Xmax, where 




VF2 _ ^2 n2 
* min ' ^ 


(27) 


and Wmin = 1-08 GeV. The values of the full integral for 
the first moment are shown in Fig. 15 as the filled (blue) 


common convention, since the quasi-elastic peak would overwhelm 
the integrals at small . 

^ The contribution from x < 0.001 is most certainly negligible for the 
higher moments. 


data points and the full systematic uncertainty due to the 
additional model uncertainty in the unmeasured region is 
indicated by the wider blue band behind the beige one. We 
also show published world data on the first moment in the 
same range. Our data are again in reasonable agreement 
with the world data (within statistical uncertainties) except 
for being slightly below the data from EGl-dvcs m as 
mentioned before; again, the difference is consistent with 
the systematic uncertainty on both experiments. At < 
0.8 GeV2, ours are the only high-precision data available so 
far, extending down to = 0.05 GeV^, where they can 
be used to test effective theories like Ghiral Perturbation 
Theory (xPT). 


We compare our results with several theoretical predic¬ 
tions and parametrizations in Fig. The black dashed- 
dotted curve indicates the extrapolation from the DIS limit 
using pQGD corrections up to third order in as, assuming 
the asymptotic value for the moment from recent publi¬ 
cations by GOMPASS [19] and HERMES [18] . We also 
show two parametrizations that connect the DIS limit with 
the real photon point. One parametrization, by Burkert et 
al. [99] (upper magenta curve), combines an estimate of 
the integral in the resonance region with a smooth function 
connecting the photon point, constrained by the Gerasimov- 
Drell-Hearn (GDH) sum rule [44] [45], with the asymptotic 
limit. The second, by Pasechnik et al. |100j (light blue line), 
includes both higher-twist terms at large and a chiral- 
like expansion at the photon point within the framework of 
an analytic perturbation theory (APT) which has been fit 
to available data, including previous (partial) results from 
EGlb [26]. Both parametrizations do a remarkably good 
job describing the world data on the first moment over the 
full range of Q^. 

We also show several predictions for the low-Q^ behavior 
of Ti on the right-hand side of Fig . [Ts] in cluding the slope 
at = 0 from the GDH sum rule [441145] (solid black line) 
and its extensions from two recent chiral perturbation theory 
calculations. The first one, by Bernard et al. [38] (narrow 
dark grey band on r.h.s.). is an expansion up to third order 
with explicit inclusion of A(1232)3/2+ isobar degrees of 
freedom. The second, by Lensky et al. [10111102] (wider 
dark green band), uses Baryon xPT including pion, nucleon 
and A(1232) degrees of freedom to calculate all moments in 
next-to-leading order (NLO). Both predictions are close to 
the GDH limit and show little sign of the observed deviation 
of the data towards less negative values as increases; 
however, they agree with our lowest three points < 0.08 
GeV^ within their statistical and systematic uncertainties. 

The higher moments Tf and Tf are also calculated in 
the same way with appropriate powers n — 3,5 (see Sec¬ 
tion 


IA I. Fig. 
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shows the results for the third moment Tf 
and the fifth moment Tf of gi from the EGlb data. These 
moments are useful for the extraction of higher twist matrix 
elements, e.g., the third moment is directly related to the 
matrix element 02 within the Operator Product Expansion. 

To calculate the extended spin polarizability 70, we in¬ 
tegrate the product of AiFi instead of gi, weighted with 
x'^. The result is multiplied by 16M'^{hc)'^a/Q^ to con- 
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FIG. 16. (Color Online) Higher moments of gf extracted 
from the EGlb data versus Q^. The third moment for the 
deuteron, Ff, is shown on the left, and the fifth moment, Fi, 
on the right. The open squares were calculated with no model 
contribution while the filled squares include model input for 
the kinematic regions with no available data. The total sys¬ 
tematic uncertainty is shown by the blue (experimental only) 
and black (experimental plus extrapolation) shaded areas. 


FIG. 17. (Color Online) Forward spin polarizability ( 70 ) 
for the deuteron versus Q^. The open squares represent the 
result using only data and the solid black circles are data plus 
model results. The shaded area close to the a;-axis is the total 
systematic uncertainty (blue for experimental only and black 
with extrapolation uncertainty included). Our model is shown 
as a red solid line. Our results are compared to three yPT 
calculations (see text) and the MAID parametrization [lO.S) 
for single pion production. 


vert to [10“^ in agreement with the definition for real 
photons. Fig. [TtI shows our result for the forward spin po¬ 
larizability 7 o for the deuteron. We compare them again to 
the xPT calculations by Lensky et al. [10111102] (upper yel¬ 
low band) and by Bernard et al. [38] (lower light blue band) 
as well as an evaluation of single pion production data by 
the MAID collaboration [103] . The xPT calculations do not 
quite reproduce the trend of the data at low Q^. 


F. Neutron spin structure functions 


Although many data sets exist for spin structure func¬ 
tions of the (bound) neutron in the deep inelastic (DIS) 
region, no un-integrated results have been published in the 
region W < 2 GeV of the nucleon resonances. This is due 
to the difficulty of reliably extracting neutron information 
from measurements that have to use nuclear targets, as ex¬ 
plained in Section []TG| As discussed in that section, we have 
attempted, for the first time, to combine our deuteron data 
with our proton fit (Section VD) and an impulse approx¬ 
imation folding prescription to access information on the 
neutron in a model-dependent way, see Figs, [island [T^ 


Our method relies on the folding prescription by Kahn 
et al. m which describes deuteron structure functions in 
terms of those of the proton and the neutron. We used this 
prescription in our fit for the asymmetries A" and A 2 for the 
neutron as described in Section [V D| In particular, for any 
set of fit parameters, we calculate both g ",52 ^nd 
combine them (following Ref. |27j) and compare directly to 


the measured gf. The parameters are optimized until the 
best possible agreement (smallest x^) is achieved. 


follows a slightly different procedure than that described in 
Ref. [27], but is similar to their “additive" method: We 
assume that any difference between the measured and the 
calculated gf is solely due to a corresponding discrepancy 
in gf at that specific kinematic point. Given that to first 
approximation 




gf^il-1.5PD)igf + gf) (28) 

we then calculate 

r,, N 71, 1 17 gf (meas) — gf (model) 

gf (meas) = gf (model) + 1 - 1 5 P^ - 

with Pd « 0.05. This method has the advantage that it 
is stable (as opposed to trying to invert the folding) and 
that it leads to a straightforward propagation of statistical 
uncertainties: 


, 717 / 7 cr(5f)(ineas) 

= 1 - l.spf 


(30) 


The systematic uncertainties are evaluat ed in th e same fash¬ 
ion as in all previous cases (see Section IV D 6), by varying 
one model input or experimental parameter in sequence and 
propagating the variation to the final result for gf (meas), 
then adding all of these variations in quadrature. The final 
results are shown in Fig. together with world data at 
higher W and both our model parametrization (red line) 
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FIG. 18. (Color online) Our resnlts for the spin structure 
function xgi of the bound nentron, extracted in the impulse 
approximation framework of Ref. |27| versus the Bjorken vari¬ 
able X for all (combined) bins (filled circles). Our model 
is shown as red lines on each plot, and the asymptotic form 
of gi(x) in the DIS region is shown as dashed blue lines. The 
shaded area at the bottom of each plot represents the system¬ 
atic uncertainty. Additional data from other experiments are 
shown as well: E154 [95] (magenta inverted triangles), HER¬ 
MES [IHIIHTI (red circles) and E142 [21] (brown triangles). 


and the DIS limit at = 10 GeV^ (blue dashed line). 
We combined our standard bins pairwise for clarity of 
presentation. 


As a next step, we can then convert the results for g" into 
the virtual photon asymmetry A", by using our models for 
F” and A^. The results are shown in Fig. 19 Overall, the 





FIG. 19. (Color Online) Ai for the bound neutron, extracted 
from our results for gi (see Fig. 181, versus W for our com¬ 
bined bins. Systematic uncertainties are shown as shaded 
area at the bottom of each plot. Our parametrized model is 
also shown as a red line on each plot. Only the data points 
with (Tstat < 0.6 and o-gys < 0.2 are plotted. The cyan dia¬ 
monds indicate data from measurements on ®He 11211881. 


agreement of the extracted results with our model is quite 
good, except at the highest where our data seem to lie 
systematically lower (a trend that can already be observed 
in the corresponding deuteron data, see Fig. 111. We direct 
the attention of the reader to the additional data points 
plotted in the last two bins (cyan diamonds); these are 
the results from the Hall-A experiment on ^He [121 US] at 
the highest attainable x in the DIS region. These data 
are consistent with our own, but with significantly smaller 
statistical uncertainties. However, no such data have been 
published for any of the lower bins. 

As a final step, we once again form various moments 
of the neutron spin structure functions (see Figs. 20 and 


211. While the advantage of using deuterium as a proxy 


for the neutron (namely, its much smaller average nucleon 
momenta and therefore less severe kinematic smearing) is 
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FIG. 20. (Color online) Fi for the neutron versus from 
data only (open blue circles) and data plus model (full blue 
circles), including the extrapolation to the unmeasured kine¬ 
matics. Also shown are phenomenological calculations from 
Pasechnik et al. | 100| (lower light blue line) and Burkert et 
at. [99] (upper magenta line), together with the xPT results 
from Lensky et al. | 1011 I102 | (wider dark green band) and 
Bernard et al. [38] (thin grey band). The GDH slope (black 
solid line) and pQCD prediction (black dotted line) are also 
shown. The right-hand side plot is a magnification of the 
low region (which is omitted from the l.h.s.). Systematic 
uncertainties of our data are shown as shaded areas at the 
bottom of the plot. Results from other experiments are also 
shown, with statistical and systematic uncertainties (added 
in quadrature) reflected in their total error bars. 


less clear in this case (since the moments integrate over 
all kinematics anyway), it is still instructive to compare our 
results to those using a ^He target as a source of polar¬ 
ized neutrons [IT]. Again, we find good agreement between 
these two experiments using different effective neutron tar¬ 
gets and with very different systematic uncertainties. We 
note that the neutron data are also well described by the two 
parametrizations [99l llOOj . while they approach the GDH 
limit above (but marginally compatible with) the Chiral Per¬ 
turbation calculations |38l llOlj . 

Figure shows the forward spin polarizability for the 
bound neutron from our data, again compared to data from 
the ^He experiment in Hall A [T^. The agreement at the 
lowest is excellent, and our data extend to slightly lower 
Q^. Once again, they show a general agreement with the 
order of magnitude predicted by x^T while exhibiting a 
distinctly different shape with Q^. 
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FIG. 21. (Color Online) Forward spin polarizability 70 for 
the neutron versus . The open squares represent the re¬ 
sult using only data and the solid black circles are data plus 
model results. The shaded area close to the a;-axis is the total 
systematic uncertainty (blue for experimental data only and 
black including the extrapolation). Our model is also shown 
as a red solid line. Our results are compared to three xPT 
calculations (see text) and to the ^He data from Hall A [11] . 


VI. CONCLUSION 

In summary, we present the final analysis of the most ex¬ 
tensive data set on the spin structure functions Ai and gi 
of the deuteron in the valence and resonance region. The 
data cover two orders of magnitude in squared momentum 
transfer, 0.05 < < 5 GeV 2 , connecting the region of 

hadronic degrees of freedom and effective theories like xT’T 
near the photon point with the regime where pQCD is ap¬ 
plicable. Our data give more detailed insight in the inclusive 
response of the deuteron in the resonance region and how 
it connects with the DIS limit. They can constrain NLO fits 
(including higher twist corrections) of spin structure func¬ 
tions extracting polarized PDFs, and they shed new light on 
the valence quark structure of the nucleon at large a;.They 
can be used to study quark-hadron duality and to extract 
matrix elements in the framework of the Operator Product 
Expansion. To facilitate such analyses, we are providing 
the raw data (with minimum theoretical bias) through the 
CLAS experimental database m as well as Supplemental 
Material for this paper [82] . 

We use our data on the deuteron, together with a detailed 
fit of the corresponding proton data, to extract bound neu¬ 
tron spin structure functions, using a convolution model and 
ignoring FSI and other binding effects. These results give 
information, for the first time, on inclusive neutron spin 
structure in the resonance region W < 2 GeV. They can 
also be used to cross check the results from ^He targets at 
high X. We find general agreement between the data from 
these rather different approaches, within the relatively larger 
statistical uncertainties of our data set. On the other hand. 
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our data cover a larger range in and W. 

Our data allow precise determinations of moments of gf 
(and g") as a function of Q^, which can be used to test 
the approach to the GDH sum rule limit, xPT phe¬ 
nomenological models, and to extract matrix elements in 
the framework of the Operator Product Expansion. We 
find that xPT describes our results for Fi only up to very 
moderate « 0.08 GeV^ (within our statistical and sys¬ 
tematic uncertainties), while there is only rough agreement 
in magnitude between xPP our data for the forward 
spin polarizability Jq. Finally, we would like to refer the 
reader to a recent analysis of the world data |104| . includ¬ 
ing the data presented here, to study the first moment of 
the difference g\ — g" and its Q^-dependence to extract 
Operator Product Expansion matrix elements. 

Further data will come from the analysis of the EG4 ex¬ 
periment with GLAS, which will extend the kinematic cov¬ 
erage of the present data set to even lower for a more 
rigorous test of xPT. Additional information on the struc¬ 
ture functions g 2 and A 2 is forthcoming once experiment 
“SANE" in Hall G and experiment “g2p" in Hall A have 
been analyzed. Finally, a complete mapping of spin struc¬ 
ture functions in the valence quark region, out to the highest 
possible X, is one of the cornerstones of the program with 
the energy-upgraded 12 GeV accelerator at Jefferson Lab. 
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